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We propose a scheme for generation of arbitrary coherent superposition of vortex states in Bose- 
Einstein condensates (BEC) using the orbital angular momentum (OAM) states of light. We devise a 
scheme to generate coherent superpositions of two counter-rotating OAM states of light using known 
experimental techniques. We show that a specially designed Raman scheme allows transfer of the 
optical vortex superposition state onto an initially non-rotating BEC. This creates an arbitrary and 
coherent superposition of a vortex and anti- vortex pair in the BEC. The ideas presented here could 
be extended to generate entangled vortex states, design memories for the OAM states of light, and 
perform other quantum information tasks. Applications to inertial sensing are also discussed. 

PACS numbers: 



Generation and manipulation of macroscopic superpo- 
sitions as well as entangled states is of paramount interest 
to the field of quantum information [ll. In this regard, 
Bose-Einstein condensates (BEC) 0, H come across 
as ideal candidates. BECs correspond to highly-coherent 
macroscopic ground-states of the confining potentials. 
Moreover, vortex states of BECs, which are topological 
states with special phase structure, have been realized ex- 
perimentally [a| ■ Stirring a BEC cloud with laser beams 
leads to the nucleation of vortex lattices in the BEC. 
These vortex states are fairly stable and could be can- 
didates for qubits in quantum information processors, if 
appropriate means to manipulate them are developed. 

In an entirely different area of optical physics jtremen- 
dous progress has been made in creation 0, H, @j 
mami julation \u%, detection [Til HI IT^. and applica- 
tion of the orbital angular momentum (OAM) states 
of light. The OAM states have a corkscrew type helical 
phase structure. To illustrate, an OAM state with an- 
gular momentum M has \£\ azimuthal phase singularities 
across a cut taken in the beam path. The sign of £ corre- 
sponds to the sense of rotation of the phase fronts around 
the beam axis. Each photon in the OAM beam carries an 
orbital angular momentum of Ti£. The quantum nature 
of these OAM states has been demonstrated recently by 
showing that a photon pair created in parametric down- 
conversion process is entangled in the orbital angular mo- 
mentum space along with the usual polarization entan- 
glement [lq . 

Excitation of vortices in BECs, using the optical vor- 
tex beams, has been proposed recently In this 
letter we introduce a scheme for creation of macroscopic 
superpositions of BEC vortex states through transfer of 
angular momentum of light from specially prepared OAM 



state superpositions. 

Generation of the superposition of gaussian beams 
with OAM states of light has been demonstrated [lflj . 
Our interest, however, lies in creating an arbitrary su- 
perposition of two counter-rotating optical vortices. 

The OAM states of light have unique amplitude and 
phase structures. To illustrate, monochromatic OAM 
beams have an azimuthal phase dependence of the type 
exp(i£(j>). Laguerre-Gaussian (LG) laser modes are an 
example of such OAM states [20|. The normalized LG 
mode at the beam waist (z = 0) and beam size wo at the 
waist is given in cylindrical coordinates (p, <f>, z) by 
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exp(-p 2 /wl)exp(i£<f>) C 1 ) 
where L l p (p) are the associated Laguerre polynomials, 

rWl \ W 1\m (1^1+ PV- rn f \ 

4>(p) = ^(-i) {p _ m)m + m)[m ,p » (2) 

wo is the beam width, p is the number of non-axial ra- 
dial nodes of the mode, and the index £, referred to as the 
winding number, which describes the helical structure of 
the wave front around a phase dislocation. For further 
discussion we consider only pure LG modes with charge £ 
and p = 0, we denote such a state of the light field by \£) 
such that (r\£) = LGg(p, <f>). Thus it can be easily seen 
that the states | +£) and | —£) , with £ being a whole num- 
ber, differ only in the sense of the winding of the phase 
either clockwise or counter-clockwise. Our aim is to cre- 
ate a general superposition of the OAM states of light 
of the kind: (a+ \£) + a_ \-£)), with |a + | 2 + |a_| 2 = 1. 
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It is well known that creation of superposition of OAM 
states of the kind J2 e eg \£) is a fairly straightforwardpro- 
cedure by using computer generated holographs or 
phase plates 0. Moreover, a sorter of these OAM states 
has also been demonstrated that can distinguish and 
separate different OAM components. Thus, by using a 
mixed OAM-state generator and a OAM-sorter in con- 
junction one can easily obtain a pure OAM state \£). 

We note that dove prisms can be used to change the 
handedness of light beams passing through them [2*lj . 
Consequently, the sense of the phase winding of an LG 
beam would be reversed as it passes through a dove 
prism. Using this, we devise a Mach-Zendcr type con- 
figuration as shown in Fig. ^ to generate a general su- 
perposition {t \£) + r \—£)) at one of the output ports of 
the interferometer. The first beam splitter is taken to 
be a special beam splitter with the ratio of i t : r for the 
transmitted and the reflected amplitudes at its output 
ports. The second beam splitter is a usual 50:50 beam 
splitter. The operation of the Mach-Zender configuration 
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FIG. 1: Scheme for creation of superposition of the OAM 
states. The first beam splitter (BS) at the input is a \f \ 2 /\i\ 2 
beam splitter as f and it are the reflection and transmis- 
sion amplitudes. The second beam-splitter is taken to be 
50:50, with the mirrors (M) perfectly reflecting. The dove 
prism (DP) transforms a right handed coordinate-system into 
a left handed one; thus it performs the operation \£) — > |— £). 
Choosing (f> — n, and discarding the output at the output port 
2, one obtains (t \+£) + r at output port 1. 

of Fig.QJcan be described throu gh t he matrix representa- 
tion of beam splitter operation |22( , such that the initial 
state (\£) , 0) T transforms into 



J_ / t\£)+f\-£) 
V2 \i(r\-t)-t\£)) 
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at the exit ports 1 and 2 of the Mach Zender Interfer- 
ometer, respectively with the choice of </> = tt. Thus, by 
ignoring port 2 and renormalizing the state from port 1 
we obtain the required superposition state t \£) + r \ —£), 
as |r| 2 + \i\ 2 = 1. In the following we present our scheme 
for transfering this optical vortex state superposition to 
BEC vortex superpositions. 

Highly detuned optical fields in a Raman configura- 
tion have been used to coherently manipulate and cre- 
ate various superpositions of different atomic levels |23| . 
Similar optical manipulation techniques exist to couple 



BEC clouds in different internal states. Moreover, the 
OAM states of light have also been shown to be useful 
for excitation of vortices in BEC through time-dependent 
linearily varying two-photon detuning [l7l | and through 
the STIRAP |2J| type scheme [l8|- In the following we 
discuss a Raman type scheme to generate superposition 
of vortex states in BEC. 

The level scheme for our model is depicted in Fig. [21 
Initially nonrotating state |0) is coupled optically via 
two Raman type configurations of the external fields, 
(Q + ,Q C ) and (r2_,f2 c ), through internal states \i) and 
\i'). The polarizations of the optical fields are taken as 
shown in the figure; thus the internal quantum numbers 
of the final states |+) and |— ) are the same. The Rabi 
frequencies f2+ and f2_ are due to the coupling of the 
BEC cloud with the two counter-rotating components of 
the special optical vortex state generated through the 
configuration discussed in Fig. ^ 




FIG. 2: The level scheme for generation of vortex state su- 
perposition. A non-rotating state |0) is coupled to the vortex 
states |+) and |— ) through the optical vortex field components 
providing coupling strengths 57+, fl- and a strong drive field 
with Rabi frequency Q c . The optical vortex beam is o"+ polar- 
ized, whereas the drive field is cr_ polarized such the hyperfine 
quantum number tuf is the same for the components |0), |±). 
The vorticity of the |±) states is ±£, respectively 



Noting that the optical fields are highly detuned, the 
intermediate states (|i) and are sparingly popu- 

lated, Thus, they can be adiabatically eliminated from 
the equations. Within this adiabatic approximation we 
can write a modified set of Ginzburg-Pitaevskii-Gross 
(GPG) [23 equations for the multi-component BEC 
trapped by the cigar-shaped trapping potential f = 
(m/2) [u>j_r + ui z z 2 ), where uj± and u> z are the trans- 
verse and longitudinal trapping frequencies respectively. 
Also to is the mass of the individual atoms in the BEC 
cloud. We note that r = \J x 2 + y 2 is the transverse ra- 
dial coordinate and <fi would be the azimuthal angle in 
the x-y plane that would be required later to describe 
the phase structure of the rotating BECs. The configu- 
ration space representations of the states |0) , |+) and |— ) 
shown in Fig. are taken to be 4 , o, v I'+ and ^ _ respec- 
tively. Thus, including the Raman- type optical couplings 
arising from the adiabatic elimination of the intermedi- 
ate state, we arrive at the modified GPG equations for 
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the three relevant components of the BEC cloud 

i h* = + ^ ( l^l 2 * + £ n * n ce- iAt ^i j 

\i=± i=± / 

i fi# ± = ± + A (|r> c | 2 * ± + fi ± f]*e lAt * ) (4) 

where = ^ + f- f i + ri{\^Q\ 2 + \^ + \ 2 + \^-\ 2 ) is the 
self-energy operator which includes the kinetic po- 
tential ('V), and the interaction energy operators for the 
BEC states. Here rj = 4irfra sc N/m signifies the strength 
of the inter-particle interactions of the iV-particle BEC 
cloud, through the scattering cross-section a sc . Using 
the LG beam mode function Eq. (JJ, the Rabi frequen- 
cies corresponding to the coupling between the optical 
vortex beam superposition a+ \£) + a_ |— £) and the BEC 
cloud can be written as 

n±(r) = a ± VL a e^ /w \^2r /w) W e ±u ^e ikz (5) 

with f2o is the usual atom-field interaction Rabi fre- 
quency. The amplitudes a± are respectively t and r 
corresponding to the first beam splitter transmission 
and reflection amplitudes of the configuration shown in 
Fig. ^ With the size of the condensate chosen to be 
much smaller than the Laguerre Gaussian beam waist, 
the exponential dependence on the radial coordinate in 
Eq. (J5J can be ignored. Now we make an ansatz that the 
topological structure of the states |0), |±) are given by 

*o(r, t) = a(t) exp[i (n/fi - n)t] i/i g (r) , 

*±(r,t) =0±(t)exp[i(d + n/%-K)t]ip v ±(r), (6a) 

where the nonrotating component tpg( r ) an d the rotating 
vortex components ip v ± (r) are given by 

V s (r) = exp{-(l/2)[(r/L ± ) 2 + (z/ L z ) 2 }} /n^L^J 2 
1 p v± (v) = (x±iyf\/^.L^ g (r). (6b) 

Here L± and L z are the size parameters of the condensate 
in the x-y plane and the z directions respectively and S 
is the two-photon detuning as shown in Fig. El Thus, the 
time dependence of the populations of different compo- 
nents (|a(t)| 2 , \(3±(t)\ 2 ) can be studied by projecting the 
rate equations (@J on to the topological states © ■ So far 
the equations are very general and no restriction exists 
on the OAM quantum number I. Hereafter, for conve- 
nience, we resort to a particular value of £ = 2. However, 
one may note that the general idea would remain valid 
for any given I. Thus, taking the projections onto the 
specified rotating or non-rotating states we arrive at 

id(t) = 3K|a(i)| 2 a(i) + oj± (a* + /3 + (t) + a*/3_(t)) 

i$±(t) = (£ + 2w x + ~ ^ lft(*)| 2 )At(*) + wl «± «(*) 
1 i=± 

(7) 



Here, the interparticle interaction strength appeares 
through the parameter k = itha sc N /[m(2ir)' i / 2 L 2 L L z \. 
This set of equations can be solved numerically, using 
the experimental parameters for a 87 Rb BEC [26j (i.e., 
lui_ = 132 Hz, a sc — 5 nm, L± = 2.35 //m, L z = 1.4 /xm), 
so that k = 422 Hz. The results of our numerical studies 
are summarized in Figures [3] and 0] We define a transfer 
function f(t) = \a(t)\ 2 - \(3+{t)\ 2 - |/3_(t)| 2 , which signi- 
fies the amount of population transferred from state |0) 
to states |±). Initially f(t = 0) = 1 as all the popula- 
tion resides in the non-rotating ground state. If complete 
transfer is achieved to an appropriate superposition of 
the |±) states then f(t) — > — 1. This transfer function is 
plotted for various values of the two-photon detuning 5 
in Fig. 01 The result to be noted is that only a continuous 
time variation of the detuning leads to complete popula- 
tion transfer at the steady state. This can be understood 
as the interaction terms in the dynamical equations effec- 
tively lead to a time-dependent stark shift of the different 
states, and only a time varying detuning maintains effec- 
tive two-photon resonance to give complete population 
transfer. It can be shown that the final state obtained 
at steady state is |+) + a* |— ) where we started with 
the optical vortex state a+ \£) + a_ \ —£), i.e., the phase 
relation between the BEC vortex components is closely 
related to the optical vortex components. In figure 01 we 
show generation of various superposition of the vortex 
states |±) for the time varying detuning. 
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FIG. 3: Transfer function f(t) for various detunings. (a) 
5 = 0, (b) 6 = 900 , (c) 8 = 380, (d) linearly varying detun- 
ing 5{t) = 3000 - 400 2 t. All detunings are given in Hz. The 
complete population transfer to the vortex state, correspond- 
ing to fit) = — 1 for sufficiently large t, is possible only with 
time-dependent detuning. 

These vortex states could be detected, for example, by 
following the proposal of Bolda and Walls 27], which in- 
volves observing a distinctive interference pattern of the 
vortex state with a non-rotating BEC cloud. However, 
such a measurement can not distinguish the sense of ro- 
tation the BEC cloud. We propose scattering of light 
incident on the cloud perpendicular to the rotation axis 
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FIG. 4: Various superpositions of the vortex states. The steady state population ratios are give in the notation \(3+\ 2 : |/3-| 2 
The vertical axis on all the plots correspond to the populations of various states and the horizontal axis is the time measured 
in seconds. 



to observe a Doppler shift of the scattered light. Such 
a measurement scheme would detect the sense of rota- 
tion of the vortex state and would collapse the state to 
either of the two counter-rotating components. We con- 
jecture that such a coherent superposition of two counter- 
rotating vortex states would be ideally suited for gy- 
ropscopy and other forms of inertial sensing. An entan- 
gled state of these vortices would also be interesting ob- 
jects to study. Further discussions of these and related 
issues will be presented elsewhere. 

To conclude, we have devised a scheme to generate su- 
perposition of two counter-rotating optical vortices. We 
have designed a Raman type scheme to transfer the op- 
tical vortex superposition onto a superposition of vor- 
tices in BEC. These macroscopic superpositions of two 
counter-rotating vortices are very interesting objects. We 
suggest a detection scheme based on rotational Doppler 
shift of a light field to sense the sense of rotation. 
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